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Abstract 

Using the gauge/gravity duality, we compute the leading order hadronic (HLO) contribution to 
the anomalous magnetic moment of muon, a^ L °. Holographic renormalization is used to obtain 
a finite vacuum polarization. We find a^ LO = 470.5 x 10 -10 in AdS/QCD with two light flavors, 
which is compared with the currently revised BABAR data estimated from e + e~ — * ir + ir~ events, 
a H LO [7T7r] = (514.1 ±3.8) x lO" 10 . 

Keywords: gauge/gravity duality, anomalous magnetic moment, muon 



1 



INTRODUCTION 



A recent high-precision measurement of 7r + 7r~ cross section by the BABAR experiment l| 
has found that the previous e + e~-based evaluation of the leading hadronic (HLO) contribu- 
tions to muon g — 2 was fairly underestimated. With this new measurement the discrepancy 
between the e + e~ and r-based results for the dominant two-pion mode is narrowed to 1.5 a 
level in the dispersion integral [2j. It is therefore quite necessary to find other means to test 
this new experimental data. 

To confirm the new BABAR data one needs to calculate it directly from QCD, for which 
Lattice QCD would certainly be the best tool. However, the state-of-the-art calculations of 
lattice QCD still suffer from a large systematic error, compared to the experimental value, 
due to difficulties in simulating light quarks sj]. 

There have been a few other attempts to calculate the HLO contribution from QCD, using 
chiral perturbation or 1/N C expansion, or from models of QCD such as the extended Nambu- 
Jona-Lasinio models m or the quark resonance model [51. In this paper we use recently 

npTiQ 

proposed models of QCD, called holographic QCD [6|, |7|, |8|, |9j, based on the gauge/gravity 
duality, to calculate the hadronic vacuum polarization contribution. The holographic models 
of QCD were found to be quite successful to account for the properties of hadrons at low 
energy and give relations to their couplings and also new sum rules [h]]. Similar attempt 



was made recently to estimate success: 
anomalous magnetic moment of muon 



Til 



y the hadronic light-by-light contribution to the 



GENERAL FORMULAS AND ADS/QCD 

Hadrons (or quarks) contribute to the anomalous magnetic moment of muon only through 
quantum fluctuations. The leading contribution is therefore the vacuum polarization correc- 
tion, shown in Fig. (TJ to the internal photon line, U^(Q 2 ), which is defined in QCD by the 
electromagnetic currents of quarks, J^ m = <?/7 M Q em <?/ (/ = u,d, s, ■ ■ • ) with the electric 
charge operator Q era , as 

i J *xe*{*\TJ^x)rjiW) = OAT - A^n^f-g 2 ) , (1) 

where rf v = diag(l, — 1, — 1, — 1). Since the (bare) vacuum polarization function is a di- 
vergent quantity, one should regularize it with a regulator and renormalize it by adding a 
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counterterm. Furthermore, if one uses the physical charge in the calculation of muon g — 2, 
the renormalized vacuum function should be modified as 

nffCQ 2 ) = nS(QV nfftf)) (2) 

so that the pole residue of the photon propagator does not get shifted. 




FIG. 1: Leading hadronic- vacuum-polarization contribution to muon g — 2. 



After the vacuum polarization function is obtained, we can calculate the hadronic leading 
order correction to the anomalous magnetic moment of muon, using the formula derived 



in 
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where Q 2 is the Euclidean momentum-squared and the kernel is given as 
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It is straightforward in holographic QCD to calculate the correlation functions of QCD 
flavor currents. According to the gauge/gravity duality, the classical gravity action of holo- 
graphic dual becomes the generating functional for the one-particle irreducible (1PI) func- 
tions of the boundary gauge theory in the large N c limit, and the bulk fields, evaluated at 
the UV boundary, serve as sources for the currents. Though the exact duality has not been 



1 An alternative form of this formula, which exactly gives the same amplitude, is presented in Ref. [T3 | 
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established for QCD yet, as a simple model, we consider a holographic model of QCD, pro- 
posed by Erlich et al. {f| for the calculation of the muon anomalous magnetic moment. Our 
analysis, however, can be easily applied to other holographic models such as Sakai-Sugimoto 
model [91] 2 

The AdS/QCD action {(J, relevant to our present computation, is given in a slice of five 
dimensional anti de Sitter (AdS) space of unit radius as 

1 

'M 

where g = det[# M jv], 9mn = (l/z) 2 r] MN with r] MN = diag(l, -1, -1, -1, -1), F R)L = 
d M R(L) N — d N R(L) M — i[R(L) M , R(L) n] with the gauge fields L M ,R M of the V(Nf) L x 
U(^V/)_r gauged-flavor symmetry in the 5D bulk and X denotes the bulk scalar fields, trans- 
forming under the gauge symmetry as X — > UlXU r . One could add higher dimensional 
operators in the AdS/QCD action ([5]) such as \X\ 2 F 2 or F 4 , but they are suppressed by a' 
and thus negligible at low energy. The bulk scalar gets a vacuum expectation value to break 
the bulk gauge symmetry down to its vectorial subgroup XJ(Nf), 

(X) = ±Mz + ^z 3 (6) 

where M = diag(m u , m^, m s , ■ ■ ■ ) is the current quark mass matrix and £ is the chiral 
condensate. 

Introducing the vector and axial vector gauge fields, V = (L — R)/2 and A = (L + R)/2, 
respectively, and working in V z = A z = gauge, we can derive the equations of motion for 
and A^ and solve them by imposing the boundary conditions for and A^ as well as for the 
bulk scalar field X, following the familiar AdS/CFT dictionary. One should note, however, 
that as far as the large- N c contributions (equivalently the meson exchanges at tree-level) to 
the vector-current correlator Uy(Q 2 ) are concerned, we are allowed to focus only on (infinite 
tower of) neutral vector mesons coupled to the vector currents as illustrated in Fig. [2j The 
values of vector meson decay constants and their masses are calculable in holographic QCD. 

Since the electromagnetic charge Q em = |T + T 3 in the case of two flavors, where the 
generators of flavor symmetry T = |l2x2 and T 3 = diag(|, — |) , the quark electromagnetic 



2 The leading hadronic contribution to muon g — 2 has been calculated in a D3-D7 setup 15j, where, 
however, scalar quarks are not decoupled. 
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V»'°'=P» to, ... <N ' 

FIG. 2: A diagram illustrating neutral vector meson exchanges giving dominant contributions to 
liy at the large N c limit. 



currents are given as 



Jem — qJq + J3 • (7) 



In the limit of exact flavor symmetry the vector-current correlator is given as 

i J d 4 xe^(0\TJ^(x)4%0)\0) = 5 a \q 2 ^ - q^)U v (-q 2 ) , (8) 
with a, b being indices of the U(2) flavor symmetry. We have then 

nfr 2 (-9 2 ) = yn y (- 9 2 ). (9) 

For Nf = 3 we have J e m~ = J3 + ^75^8 an d n e m~ (— q 2 ) = |lly(— q 2 ) in the exact limit of 
the U(3) flavor symmetry. 

HOLOGRAPHIC VACUUM POLARIZATION 

Since we are interested in two-point functions of flavor currents, we need to keep terms 
quadratic in the bulk gauge fields for a given vacuum solution for (X) in the action (JHJ). If 
we neglect the higher order terms like |X| 2 F 2 , the vector gauge fields do not couple to the 
chiral condensate S but may couple to the current quark mass operator, M. However, at 
the leading order, the neutral vectors do not couple to X at all and have same mass among 
them, since (X)] = 0. 

The equations for vectors are given at the leading order as 

(d z - z d z + ^)vi ll (q,z) = 0, (10) 

where V± (q, z) = v°(q)V±(q, z) is the transverse component of the Fourier transform of the 
vector field V^(x,z) = V°(x, z)T a . Imposing the boundary condition d z V^(q, z m ) = and 
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V®(q, e) = v*(q), one gets the solution to Eq. flTUJ) given by Bessel functions. Putting the 
solution into the action, we are left with the boundary action 



d A q 



4 A 4 



1 V±(q,z)d z V ± (q,z) 

gl z 



v^{q) 



(11) 



which becomes the generating functional for the 1PI correlation functions of the QCD flavor 
currents, [v^] = S^v^] . 

From Eq. (jSJ) we immediately read off the transverse part of the vector current two-point 
functions, 



S 2 S eS [v lJ 



% I d 4 xe tqx {0\TJ^(x)4 u (0)\0) 



(12) 



5v^(-q)5v vb (q) 

Using Eq. (IT2l and the action (1111) . we find the vector current correlator Uy(q 2 ) in terms of 
the five dimensional gravity dual, 



gl 



d z V±(q,z) 
q 2 z 



1 1 J (qz m )Y (qe) - Y (qz m )J (qe) 



(13) 



9l <?e Jo{qz m )Yi(qe) - Y (qz m )Ji(qe) ' 

Note that this Uy (q 2 ) behaves like In e as the UV cutoff e — * 0, so it needs to be renormalized. 
To make this point clearer, we may expand Eq. f|T3|) perturbatively in power of e to obtain 
the asymptotic form of Uv(q 2 ), 



iW) = - 

#5 



7T Y (qz r . 
2 Jo(qz r , 



+ 7-ln2 + lnge + C(e 2 



(14) 



where 7 is the Euler constant. In order to renormalize the vacuum polarization function, we 



, |l7| , which is appropriate in 



shall adopt a procedure of the holographic renormalization [if], 
our holographic approach. The holographic renormalization can be performed by introducing 
a counterterm at the UV boundary, respecting all the residual 4d-gauge symmetry, which 
subtracts only the divergent term in the vacuum polarization (fl4l . just like the familiar 
minimal subtraction scheme. We find the counterterm appropriate to renormalize Uy(q 2 ), 



d 4 x lne/i^j tr [F^ 



where /1 is a renormalization scale and = d^Vu — dyV^ — i\V^V v 
contributes to the vacuum polarization function 
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This counterterm 
(16) 



Adding this counterterm contribution to Eq. ()14p . we obtain the renormalized expression 
for the vacuum polarization function, 



n r y en ( 



lim[n y (g 2 )+n c y (g 2 )] = ^ 



nY {qz m ) . . g 

■- — -+7- In 2 + In - 

2 J {qz m ) fx 
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We may derive another renormalized expression for Ily alternative to Eq. (|T71) by expand- 
ing Eq. ([HP in terms of an infinite tower of vector meson exchanges. This is possible due to 
the existence of poles in time-like momentum region provided by zeros of the denominator 
of ily(g 2 ) of Eq. ({TBI as Jo(qz m ) = 0. Let us first consider a normalizable solution ?/y n 
to Eq. ffTUj) with q 2 replaced by M Vn , satisfying the boundary condition, d z ipy n (z m ) = 0, 
ipv n { € ) = 0) normalized by j Zm — {i>v n { z )) 2 — 1- The vector current correlator ITy(g 2 ) can 
then be expressed in terms of the infinite set of ipy n along with their poles to be with e — > 

_[<k(e)/e] 2 



oo 

n v(? 2 ) = 4£ 



(18) 



- (q 2 - Ml)Ml ' 

where the dot denotes the derivative with respect to z, ip = d z ip. From this, we read off 
the decay constants of the vector mesons, F Vn , defined by (0| J{T(0)|V^) = F Vn 5 ab e^ with the 
polarization vector e M , 

TW„(e) 




(19) 



<E^0 



One can easily show that, in the expression of Eq. (fT8"|) . the UV divergence of Eq. (fill) is 
now embedded in IIy(0) as 



n v (o) 



F 2 



oo 

5^ Mi 

n=l Vn 



- , Z m 
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(20) 



Decomposing Tl v (q 2 ) as IIy(g 2 ) = ITy(O) + [Ily(g 2 ) — IIy(0)] and replacing the first IIy(0) 



with the logarithmic divergent term of Eq. ff20|) . adding the counterterm contribution of 
Eq. (|T6|) to this ily(g 2 ), we arrive at another form of the renormalized ITy(g 2 ), Uy n (q 2 ), 



1 a F 

IiT{q 2 ) = lim[ny(g 2 ) + U v (q 2 )] = - - hi z m/ i + ^ ^ -jLJ 

in which (Fv n /M v )\ e =o is finite: 



71=1 v 



(21^ 



e=0 
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M 2 

Vn 



V2 



t=0 



where we have used Eq. ( 1191 and 



Vv„(z) 



£f 2 Z m My n Ji(2 m My n ; 
V^Jl^MyJ 



(22) 
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HOLOGRAPHIC ESTIMATION OF a 



In this section we calculate the leading order hadronic (HLO) contribution to the muon 
anomalous magnetic moment a^ LO in holographic QCD. 

We should first note that the holographic description of the present model breaks down at 
a high-energy scale above which the underlying stringy ingredients would not be negligible. 
Recall, furthermore, that the vector current correlator Uy includes the UV divergence which 
has been renormalized and converted into the renormalization scale based on the holographic 
renormalization scheme. This implies that the HLO calculation of Uy is potentially sensitive 
to the UV scale/renormalization scheme. From these points of view, it seems reasonable to 
truncate the infinite tower of vector mesons at a finite level n. The simplest way to do is 
to use Eq. fT2~Tl) and expand it in powers of (q 2 /M Vn ), and keeping only a few excited states. 
Examining the effects from higher resonances on the integral kernel of Eq. ([3]) [See Fig. [3] , 
we find that truncating the tower at n = 4 suffices within about 1% accuracy to saturate 
the full contributions. We therefore truncate at the 4-th (n = 4) excited states and derive 
an expression of Uy n (q 2 ) valid for a low-energy region My < q 2 <C {M v& , ■ • • }, to get 

nrV) = In W + £ {q2 ful)M v + °foV«» • (24) 
In order to utilize the formula Eq. ([3]), we subtract Hy n (0) to obtain 

n v (Q 2 ) = nng 2 ) - n r y cn (o) , 

= ± wSjw +om)) ' (25) 

where Q 2 = —q 2 is Euclidean momentum-squared. Note that the dependence of the renor- 
malization scale \x has been removed. 

Let us now evaluate fly of Eq. (125]) numerically. Looking at Eq. (1251) . to this aim, we 
see that all we need to do is calculate the values of the vector meson poles My n and their 
pole residues (Fy n /M Vn ) for each n = 1,2,3,4. We first notice that the mass hierarchy of 
vector mesons are completely determined by locations of zeros given as Jo(z m My n ) = as 
seen from Eq. (TIB"]) . Especially, for n — 1, 2, 3, 4, we have 

z m M Vl = 2.405 , z m My 2 = 5.520 , 

z m M V3 = 8.654 , z m M V4 = 11.792 . (26) 
8 
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FIG. 3: Comparison of the integral kernel / (Q 2 )f[ e m{Q 2 ) which has a peak around Q 2 = m 2 ^. The 
dashed curve corresponds to the result including full contributions from the infinite tower of vector 
mesons, while four bold curves are obtained by integrating out the infinite tower at the levels of 
n = 1,2,3,4. The dashed curve is almost (within about 1% deviation) reproduced when n = 4. In 
the plot the number of flavors Nf is taken to be 2. 




Therefore, once the lowest pole is identified with the p meson pole, i.e., My 1 = M p = 

ri 

0.775 [GeV] [18| which fixes z m = 3.101 [GeV] -1 , we can obtain a set of the vector meson 
masses (in unit of GeV); 

M Vl = 0.775 (fit) , M V2 = 1.779 , 

My 3 = 2.788 , M V4 = 3.802 , (27) 



in which the values of My 2 ,My 3 are compared to the empirical ones [18j], M P2 \ exp = 1.45 
and M P3 | cxp = 1.70 in unit of GeV. We next turn to calculating the pole residues F Vn /M Vn 
expressed as in Eq. ( I2"2"j) . Note that the overall factor of Eq. (1221) is determined by 

considering a fit with the expression derived from the operator product expansion of Uy(Q 2 ), 
and it is fixed as 

H = v&- (28) 

For n = 1,2,3,4, the values of Fy n /M v of Eq. (122]) are then calculated up to the overall 
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factor as 



||- ■ ( v/iVcA^T 1 = 1-133 , || • ( v/Ayi^)- 1 = -0.753 

/ /nr /10 _2\-1 _ n .no F Vi 



3 ■ ( v/iVc/12^)" 1 = 0.602 , -I- ■ ( v^yi^)- 1 = -0.516 . (29) 

V3 v& 

Now we are ready to evaluate the HLO contribution to the muon anomalous magnetic 
moment a^ L °. We shall focus on a two-flavor (iVj = 2) HLO contribution, coming from the 
lightest two flavors (it, d), and then compare the predicted value of a^ LO with that estimated 
from the new BABAR data on the e + e~ — > tt + tt~ events detected in a range of center of 
mass energies 2m n < ^/s < 1.8 GeV [l|. Taking into account Eq. and substituting fly 
of Eq. ( 1251) into the formula Eq. (J3j), using the theoretical values for the vector meson poles 
and their pole residues listed in Eqs.( l271) and ( 1291 . we obtain 

< L °lSs/ 2 QCD = 470.5 xlO- 10 , (30) 

which agrees, within 30% errors, with the currently updated value [l| 

< L °MIbabar = (514.1 ± 3.8) x lO" 10 . (31) 

We expect that the 1/N C corrections together with the isospin-breaking corrections could 
make up for the discrepancy between the values of Eq. fl30|) and Eq. fl3TT) . 



DISCUSSION 



We have calculated the hadronic leading correction (HLO) to the muon magnetic anomaly 
with two light flavors. The result is compared with a recently updated HLO value [2j, 
estimated by including not only the BABAR new e + e~ — > 7r + 7r~ data [l| but also the other 
experimental data involving other inclusive decay processes in which multi-hadrons more 
than two its, such as 27r°7r + 7r~, 2tt + 2tc~, are included in the final states, as well as some 
exclusive decay processes. From a theoretical point of view, those multi-hadronic decaying 
processes are thought to be dominated by higher order contributions of meson-loops in 
the large N c expansion. However, those higher order corrections have not been achieved 
in the present holographic framework which is based on the holographic correspondence, 
established only at the large iV c limit. It is not adequate, therefore, to compare the result 
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obtained from our present holographic calculation with that deduced from data involving 
such those decaying processes. 

Similarly, we have not compared with the value estimated based on lattice calculation [lc 
with the staggered fermions of Nf = 3, because the value of Ref. [3] includes a part of the 
meson-loop contributions. 

We suspect that incorporation of 1/N C corrections as done in Ref. jl9] might make it 
possible to compare with these references. 
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